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Nuclear modification of structure functions in lepton scattering 



S. Kumano a * 

a Department of Physics, Saga University, Saga, 840-8502, Japan 

We discuss nuclear structure functions in lepton scattering including neutrino reactions. First, the determina- 
tion of nuclear parton distribution functions is explained by using the data of electron and muon deep inelastic 
scattering and those of Drell-Yan processes. Second, NuTeV sin 2 8w anomaly is discussed by focusing on nuclear 
corrections in the iron target. Third, we show that the HERMES effect, which indicates nuclear modification 
of the longitudinal-transverse structure function ratio, should exist at large x with small Q 2 in spite of recent 
experimental denials at small x. 



1. Introduction 

Modification of nuclear structure functions or 
nuclear parton distribution functions (NPDFs) 
is known especially in electron and muon scat- 
tering. In neutrino scattering, such nuclear ef- 
fects have not been seriously investigated due to 
the lack of accurate deuteron data. Nuclear ef- 
fects in the PDFs have been investigated mainly 
among hadron structure physicists. However, de- 
mands for accurate NPDFs have been growing 
from other fields in the recent years. In fact, one 
of the major purposes of this workshop P is to 
describe neutrino-nucleus cross sections for long 
baseline neutrino experiments, so that neutrino 
oscillations could be understood accurately |2i;ij . 

In the near future, neutrino cross sections 
should be understood within a few percent level 
for the oscillation studies |2I3| . Because typi- 
cal nuclear corrections in the oxygen nucleus are 
larger than this level, they should be precisely 
calculated. In low-energy scattering, the nuclear 
medium effects are discussed in connection with 
nuclear binding, Fermi motion, short-range cor- 
relations, Pauli exclusion effects, and other nu- 
clear phenomena. In this paper, the nuclear cor- 
rections are discussed in the structure functions 
and the PDFs by focusing on the high-energy re- 
gion. These studies are important not only for 
the neutrino studies but also for other applica- 
tions. For example, they are used in heavy-ion 

*kumanos@cc. saga-u.ac.jp, http://hs.phys.saga-u.ac.jp 



physics 0] for understanding accurate initial con- 
ditions of heavy nuclei, so that one could make a 
definitive statement, for example on quark-gluon 
plasma formation, in the final state. They could 
be also used in understanding nuclear shadowing 
mechanisms [3]. 

In this paper, recent studies are explained on 
the nuclear effects which are relevant to high- 
energy neutrino scattering. First, a recent NPDF 
\ 2 analysis is reported. Although the unpolar- 
ized PDFs in the nucleon have been investigated 
extensively 0, the NPDFs are not well studied. 
However, there are some studies to obtain op- 
timum NPDFs by using a simple parametriza- 
tion form and nuclear scattering data |7I8| . We 
explain the current situation. Second, NuTeV 
sin 2 9w anomaly [§] is investigated in a conserva- 
tive way, namely in terms of nuclear corrections 
|1()I11I12I13| . The NuTeV collaboration obtained 
anomalously large sin 2 9w . Before discussing any 
new physics mechanisms |14| . we should exclude 
possible nuclear physics explanations. In par- 
ticular, the used target is the iron and it may 
cause complicated nuclear medium effects. Third, 
the HERMES effect which is nuclear mod- 
ification of the longitudinal-transverse structure 
function ratio, is investigated in a simple convo- 
lution model. It is intended to show that such 
an effect should exist in the medium and large 
x regions |lf>| in spite of recent experimental de- 
nials at small x |17ll8j . In particular, the nucleon 
Fermi motion in a nucleus could play an impor- 
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tant role for the nuclear modification. 

This paper consists of the following. In sec- 
tion |21 global NPDF analysis results are shown 
The sin 2 0w anomaly topic is discussed in section 
EJ The HERMES effect is explained in section H 
The results are summarized in section 

2. Nuclear parton distribution functions 

The determination of the NPDFs is not still 
satisfactory in comparison with the one for the 
nucleoli. It is partly because enough data are 
not obtained for fixing each distribution from 
small x to large x. For example, various scal- 
ing violation data are not available, unlike the 
HERA data for the proton, at very small x for 
fixing gluon distributions. However, the determi- 
nation of nuclear PDFs has been awaited for de- 
scribing high-energy nuclear scattering phenom- 
ena, including neutrino-nucleus and heavy-ion re- 
actions. Some efforts have been made to provide 
practical parametrizations for the NPDFs, such 
as the ones by Eskola, Kolhinen, Ruuskanen, Sal- 
gado [7] and the ones by the HKM analysis 
In the following, the NPDFs are discussed based 
on the latter study in Ref. [EJ. 

First, the parametrization form should be se- 
lected. From the studies of nuclear F^ struc- 
ture function ratios F^/F^ 3 , one knows the exis- 
tence of shadowing phenomena at small x, anti- 
shadowing at x ~ 0.2, depletion at medium x, 
and then a positive nuclear modification at large 
x. In order to express such x dependence, the fol- 
lowing functions are used for the initial NPDFs 
at Q 2 =l GeV 2 : 

ff(x,Ql)=w i (x,A,Z)f i (x,Q 2 ), 

di(A, Z) + biX + c-iX 2 + diX 3 

x (l^xj* ■ w 

Here, Z is the atomic number, A is the mass num- 
ber, and the subscript i indicates a distribution 
type: i—u v , d v , q, or g. The functions ff and 
fi are the PDFs in a nucleus and the nucleon, 
respectively, so that the weight function Wi indi- 
cates nuclear medium effects. The nuclear mod- 



ification Wi — 1 is assumed to be proportional to 
1 — l/A 1 / 3 , and its x dependence is taken to be a 
cubic functional form with the 1/(1 — x)P* factor 
for describing the Fermi-motion part. The pa- 
rameters a,-, hi, Ci, and di are determined by a \ 2 
analysis of experimental data. 

Although the flavor dependence of the anti- 
quark distributions is known in the nucleon 
the details of nuclear antiquark distributions can- 
not be investigated at this stage. Therefore, 
flavor symmetric antiquark distributions are as- 
sumed in the parametrization. 

The electron and muon deep inelastic exper- 
imental data and Drell-Yan data are fitted by 
the NPDFs in Eq. Q. The initial NPDFs are, 
of course, evolved to various experimental Q 2 
points, and \ 2 values are calculated in compar- 
ison with the data for electron and muon deep 
inelastic scattering and Drell-Yan processes: 

x 2 = Y {R i z3i ] . (2) 

j W I 

Here, R is the ratio FA/FA' or ct^y/^dy- These 
structure functions and the DY cross sections are 
calculated in the leading order. The experimental 
error is given by systematic and statistical errors 
as ( a p ta ) 2 = {af s ) 2 + {af at ) 2 . The first version 
was published in 2001, and then the research is 
in progress by including the Drell-Yan data. We 
discuss the obtained NPDFs by these analyses. 

Obtained optimum distributions are shown for 
the calcium nucleus at Q 2 = l GeV 2 in Fig. ^ 
The solid, dashed, and dotted curves indicate the 
weight functions for the valence-quark, antiquark, 
and gluon distributions. The valence distribution 
is well determined in the medium x region, but it 
is difficult to determine it at small x although it 
is constrained by the baryon-number and charge 
conservations. In fact, it will be one of the NuMI 
projects [201 to determine the valence-quark (F 3 ) 
shadowing in comparison with the antiquark (F2 ) 
shadowing by neutrino-nucleus scattering. On 
the other hand, the antiquark distribution is well 
determined at small x; however, it cannot be fixed 
at medium x (x > 0.2) in spite of the momentum- 
conservation constraint. Because this is the lead- 
ing order analysis, the gluon distribution is not 
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Figure 1. Obtained weight functions for the cal- 
cium nucleus at Q 2 —l GeV 2 . 



fixed in the whole x region. 

The obtained NPDFs are available at the 
web site http://hs.phys.saga-u.ac.jp/nuclp.html, 
where computer codes are available for calculat- 
ing the distributions at given x and Q 2 for a re- 
quested nucleus. The nuclear type should be in 
the range, 2 < A < 208, in principle, because the 
analyzed nuclei are in this range. However, the 
distributions could be also calculated for larger 
nuclei (A > 208) because variations of the NPDFs 
are rather small in such a large- A region. If one 
wishes to use an analytical form, the distributions 
at Q 2 = l GeV 2 are provided in the appendix of 
Ref. 8 . After the first version was published, a 
new analysis has been investigated. The second 
version will become available within the year of 
2003. 

3. A nuclear physicist's view of sin 2 8w 
anomaly 

The NuTeV collaboration announced that 
their measurement of the weak mixing angle 
sin 2 9w is significantly different from collider 
measurements. If the neutrino-nucleus scatter- 
ing data are excluded, a global analysis indicates 



^2 non— shell 



0.2227±0.0004 21 . On the other 



hand, the NuTeV reported 

sin 2 9 w = 0.2277 ± 0.0013 (stat) ± 0.0009 (syst) , 

(3) 



by using their neutrino and antineutrino scatter- 
ing data. 

Because it is one of the important constants in 
the standard model, we should find a reason for 
the discrepancy. Although it may indicate the ex- 
istence of a new mechanism |14| , we should seek a 
conservative explanation first. In particular, the 
NuTeV target is the iron nucleus so that nuclear 
medium effects might have altered the sin 2 9w 
value |10I11I12I13| . In the following, we explain 
nuclear effects on the sin 2 9w determination. 

The neutrino and antineutrino cross section 
data are analyzed by a special Monte Carlo code, 
so that it is not theoretically straightforward to 
investigate a possible explanation. In order to 
simplify the investigation, we study nuclear ef- 
fects on the Paschos-Wolfenstein (PW) relation, 
which is considered to be "implicitly" used in the 
NuTeV analysis. The PW relation |22] is given 
by the ratio of charged current (CC) and neutral 
current (NC) cross sections: 



R~ = 



T isN 

NC 



'NC 



T vN 

'CC 



T vN 

'CC 



sin t>w ■ 



(4) 



This relation is valid for the isoscalar nucleon; 
however, corrections should be carefully investi- 
gated for the non- isoscalar iron target. 

If the relation is calculated for a nucleus in the 
leading order of a s , we obtain ^2] 



r: 



vA 
'NC 



vA 
'NC 



vA 
J CC 



vA 
7 CC 



= {1 - (1 - y) 2 } [ («i - u\){u${x) + c£(x)} 

/[d*{x) + 4(x) - (i - yf + 4{x)} ; 



(5) 



where the valence quark distributions are de- 



fined by 



The couplings are 



expressed by the weak mixing angle as ul = 
1/2- (2/3) sin 2 9w, u R = -(2/3) sin 2 6 w , d L = 
-1/2 + (1/3) sin 2 9 w , and d R = (1/3) sin 2 9 w . It 
is known that the nuclear distributions are modi- 
fied from those for the nucleon. The modification 
for and could be expressed by the weight 
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functions w Uv and Wd v at any Q 2 : 



u A (x) = w Uv {x,A,Z) 



Z u v (x) + N d v (x) 
A ' 



d - { x)= Wdv (x,AZ) Zd ^ + A Nu ^\ (6) 

although Wi in section 2 is defined at fixed Q 2 
(—Qq). Here, u v and d v are the distributions in 
the proton, and iV is the neutron number. 

In order to find possible deviation from the PW 
relation, we first define a function which is related 
to the neutron excess in a nucleus: e n (x) — [{N — 
Z)/A](u v — d v )/(u v + dg), and then a difference 
between the weight functions is defined by 



e v (x) 



w dv (x, A, Z) - w Ui] (x, A, Z) 
w d „ (x, A, Z) + w Uv (x, A, Z) 



(J) 



Furthermore, there are correction factors associ- 
ated with the strange and charm quark distri- 
butions, so that we define e s and e c by e s = 
s A /[w v (u v + d v )} and e c = c A /[w v (u v + d v )} with 

w v = (wd v + "-OA 

Neutron-excess effects are taken into account in 
the NuTeV analysis as explained by McFarland et 
al. 1 an( i they are also investigated by Kulagin 
|13| . The strange quark (e s ) contribution is small 
according to Zeller et al. 23 , and it increases 
the deviation. Here, we investigate a different 
contribution from the £ v {x) term U^J. Writing 
Eq. (J5J) in terms of the factors, e n , e v , e s , and 
e c , and then expanding the expressions by these 
small factors, we obtain 

f 



R A = - - sin"V w 



£v{x){ ( - 



2 ^i + (i- y y 

sin v w 1 



i-(i-y) 2 



-sin V w 



+ 0(e 2 v ) + 0(e n ) + 0(e s ) + 0(e c ). 



(8) 



Because only the e v contribution is discussed in 
the following, other terms are not explicitly writ- 
ten in the above equation. This equation in- 
dicates that the observed sin 2 9w in neutrino- 
nucleus scattering is effectively larger if the ratio 
is calculated without the e v correction. 

The nuclear modification difference e v (x) is not 
known at all at this stage. We try to estimate it 



theoretically by using charge and baryon-number 
conservations: Z = J dx A^2 q e q (q A — q A ) and 

A = J dxAJ2 q (l/3) (q A - q A )- These equations 
are expressed by the valence-quark distributions, 
then they becomes 



dx (u v + d v )[ Aw (x)e n (x)}=0, (9) 

dx(u v + d v ) [Aw v {I — 3 e n (x) } 

-i« o e„(a;){3-e n (a;)}]=0 ) (10) 



where Aw v is defined by Aw v — w v — 1. These 
equations suggest that there should exist a fi- 
nite distribution for e v (x) due to the charge and 
baryon-number conservations. However, there is 
no unique solution for these integral equations, so 
that the following discussions become inevitably 
model dependent. 

We provide two examples for estimating the or- 
der of magnitude of the effect on sin 2 6w- First, 
the integrands of Eqs. @ and are assumed 
to vanish by neglecting the higher-order terms 
0(e v e n ): 



case I: e v {x) = —e n (x) 



Aw v (x) 
w v (x) 



(11) 



Second, the x 2 analysis result [Sj, which is ex- 
plained in sectional could be used for the esti- 
mation: 



case 2: e v (x) 



w d Jx) - w Uv (x) 



X analysis 



(12) 



These two descriptions are numerically estimated 
and the results are shown at Q 2 =20 GeV 2 in Fig. 
121 The solid and dashes curves indicate the case 
1 and 2, respectively. 

In the first case, the function e v is directly pro- 
portional to the nuclear modification Aw v (x), so 
that it changes the sign at x ~ 0.2. In comparison 
with the NuTeV deviation 0.005, which is shown 
by the dotted line, is of the same order of 
magnitude. On the other hand, the second one 

(2) 

e v is rather small. This is partly because of the 
assumed functional form in the \ 2 analysis 
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Figure 2. The function e v (x, Q 2 ) is estimated by 
two different descriptions at Q 2 =20 GeV 2 . 



which was not intended especially to obtain the 
nuclear modification difference e v . Because the 
distributions are much different depending on the 
model, numerical estimates are merely considered 
to be an order of magnitude estimate. 




1 10 100 



Q 2 (GeV 2 ) 

Figure 3. Contributions to sin 2 0w are calculated 
by taking the x average and they are shown by the 
solid curves. The dashed curves are calculated by 
taking the NuTeV kinematics into account |12| . 



As far as we see Fig. |2J our mechanism seems to 
a promising explanation for the NuTeV anomaly. 
If a simple x average is taken for the e v contri- 
bution to the sin 2 9w determination, we obtain 
the solid curves in Fig. [3J and they are of the 



order of the NuTeV deviation. However, the sit- 
uation is not so simple. Although the function 
Sv is large in the large x region in Fig. |3 few 
NuTeV data exist in such a region. It means that 
the e v contribution to sin 2 9w could be signifi- 
cantly reduced if the NuTeV kinematics is taken 
into account. The guideline of incorporating such 
experimental kinematics is supplied in Fig. I of 
Ref. [2^1. The distribution e v could be effectively 
simulated by the NuTeV functions, — d v v and 
d v v — although physics motivation is completely 
different. Using the NuTeV functionals |23I24| . 
we obtain the dashed curves in Fig. [3] Because 
of the lack of large x data, the contributions are 
significantly reduced. 

The mechanism due to the nuclear modification 
difference between u v and d v could partially ex- 
plain the NuTeV deviation, but it is not a major 
mechanism for the deviation according to Fig. |2| 
However, the distribution e v itself is not known 
at all, so that it would be too early to exclude 
the mechanism. On the other hand, it should be 
an interesting topic to investigate e v experimen- 
tally by the NuMI [201 and neutrino-factory [2"B] 
projects. 

4. HERMES effect 

The HERMES effect indicates nuclear modi- 
fication of the longitudinal-transverse structure 
function ratio R(x,Q 2 ). It was originally re- 
ported at small x (0.01 < x < 0.03) with small 
Q 2 (0.5 < Q 2 < 1 GeV 2 ) in the HERMES paper 
in 2000 PH|- There are theoretical investigations 
on this topic in terms of shadowing 26 and an 
isoscalar meson |27) . 

This interesting nuclear effect is, however, not 
observed in the CCFR/NuTeV experiments |17| . 
Although the CCFR/NuTeV target is the iron nu- 
cleus, observed R values agree well with theoreti- 
cal calculations for the nucleon in the same kine- 
matical region with the HERMES. Furthermore, 
a more careful HERMES analysis of radiative cor- 
rections showed such modification does not exist 
anymore |18j . 

Considering these experimental results, one 
may think that such a nuclear effect does not 
exist at all. However, we point out that the ef- 
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feet should exist in a different kinematical region, 
namely at large x with small Q 2 JS]. The ex- 
istence of a nuclear effect in R(x,Q 2 ) is impor- 
tant not only for investigating nuclear structure 
in the parton model but also for many analyses of 
lepton scattering data. For example, the SLAC 
parametrization in 1990 |28j has been used as a 
popular one. However, the data contain nuclear 
ones, so that one cannot use it for nucleon scat- 
tering studies if large nuclear effects exist in the 
data. Because of the importance of R(x,Q 2 ) in 
lepton scattering analyses, we investigate the pos- 
sibility of nuclear modification theoretically. 

The structure function for the photon polariza- 
tion A is W£' N = el* e\ W^ N , so that longitu- 
dinal and transverse ones are defined by w£ 
{W*f + W*f)/2 and w£> N 
N and A denote the nucleon and a nucleus, re- 
spectively. Lepton-hadron scattering cross sec- 
tion is described by a lepton tensor multiplied by 
a hadron tensor W^ u . In the electron scattering, 
the tensors for the nucleon and a nucleus are given 
by 



Pf W£ v = W£. They are explicitly written 

as Pr = -(1/2) (<r -p% pZ/fi) and P^ = 
-p\l{2p\) (g» v - 3p^/p 2 ). Instead of W 1 
and W2 structure functions, the functions F\ and 
F 2 are usually used: Fp N = ^Jtf~Wp N and 

F 2 A ' N = {P A , N ■ q/Jpi~ N )W^ N . Then, the Ion- 
gitudinal structure function is given by 

Q 2 



F 



£' N {x A , N1 Q 2 )=[l + 



zpA,N r 
F 2 0; 



-2x AN F^ N (x A:N1 Q') 1 (15) 

where xa — Q 2 /{2pA • q) and xn — Q 2 /{2pn ■ <?)• 
The ratio R A of the longitudinal cross section to 
the transverse one is expressed by the function 
R A (x A ,Q 2 ): 



A,N 



Wq' N . Here, 



Ra{xa,Q 



, f£{x a ,Q 2 ) 



2x a F 1 a (xa,Q 2 



(16) 



Applying the projection operators and 
P£ v to Eq. (UU), we have 



Wfr N (p A>N ,q) = -W^ N (p A:N ,q) (9^ ^ 2x A F 1 A (x A ,Q 2 ) = J d i p N S(p N )z 



M N 



+ Wt' N {p A . N ,q) 



J A,N^ l J A,NV 
* A N 



(13) 



1 + ^t) 2XNF " (XN ' Q2) + % F " 



where p^ = p^ — (p ■ q)q^/q 2 - In terms of 
these structure functions, the longitudinal and 
transverse structure functions are given 
Wr N and W£> N = (1 + v\JQ 2 )Wp N - W*' N 
by taking the nucleus or nucleon rest frame. Here, 
v A = v 1 and the photon momentum in the nu- 
cleon rest frame is denoted {vn^n) with v N = 

(pn ■ q) 2 /p 2 N - 

We use a conventional convolution description 
for nuclear structure functions: 



~— r L (xn,Q 2 ) 
(17) 



F a (xa,Q 2 ) = d i p N S{ PN )z 



M. 



N 



W*(p A ,q) 



d 4 p N S(p N )W'l( PN ,q) , (14) 



where pi^ is the nucleon momentum and S(pn) 
is the spectral function which indicates the nu- 
cleon momentum distribution in a nucleus. In 
order to investigate the longitudinal and trans- 
verse components, we introduce projection op- 
erators which satisfy P\ V W A — W A and 



+- ^AfS{x n ,Q 2 ) + ?i±2x N F?{x N M 2 ) 
Pn J Pn 

(18) 

These results are interesting. The transverse 
structure function for a nucleus is described not 
only by the transverse one for the nucleon but also 
by the longitudinal one with the admixture coef- 
ficient p^j_/(2p^). The pn± is the nucleon mo- 
mentum component perpendicular to the photon 
one q. Equations (|17f) and (|18f) indicate that the 
transverse-longitudinal admixture exists because 
the nucleon momentum direction is not necessary 
along the virtual photon direction. 

These expressions are numerically estimated 
for the nitrogen nucleus by taking a simple shell 



G 



model for the spectral function with density de- 
pendent Hartree-Fock wave functions. Parton 
distribution functions are taken from the MRST- 
1998 version and the nucleonic R(x,Q 2 ) is taken 
from the SLAC analysis in 1990 [21]- The 
nitrogen- nucleon ratios Rn n / Rn are shown at 
Q 2 =l, 10, 100 GeV 2 by the solid curves in Fig. 
0] In order to clarify the admixture effects, the 
ratios are also calculated by suppressing the p^± 
terms, and the results are shown by the dashed 
curves. In addition, the nuclear modification is 
calculated at Q 2 =0.5 GeV 2 by using the GRV94 
parametrization for the PDFs. It is intended 
to find the modification magnitude at smaller 
Q 2 , where JLab experiments could possibly probe 
|29j. Because the admixture is proportional to 
Pn±/(^Pn) ~ Pn-l/Q 2 ' * ne modification effects 
are large at small Q 2 (=0.5 — 1 GeV 2 ) and they 
become small at large Q 2 . However, the modifi- 
cation does not vanish even at Q 2 =100 GeV 2 due 
to the Fermi-motion and binding effects which are 
contained implicitly in the spectral function. 




without L-T mixing 
06 08 "~ 



Figure 4. The nitrogen-nucleon ratio Rh^/Rn is 
shown at Q 2 =0.5, 1, 10, and 100 GeV 2 . The solid 
curves are the full results and the dashed ones are 
obtained by terminating the admixture effects. 



In this way, we found that the nucleon Fermi 
motion, especially the perpendicular motion to 
the virtual photon direction, and the nuclear 
binding give rise to the nuclear modification of the 
longitudinal-transverse ratio R(x,Q 2 ). However, 
nuclear modification of R in the large x region 



with small Q 2 has not been investigated exper- 
imentally. The situation is clearly illustrated in 
Fig. 3 of Ref. JJjj where the data does not ex- 
ist at x = 0.5 with Q 2 wl GeV 2 . We hope that 
future measurements, for example those of JLab 
experiments |29j . are able to provide clear infor- 
mation on the nuclear modification in this region. 

5. Summary 

Current neutrino scattering experiments are 
done with nuclear targets, so that precise nuclear 
corrections should be taken into account in order 
to investigate underlying elementary processes, 
for example neutrino oscillation phenomena. In 
this paper, the discussions are focused on high- 
energy reactions. 

First, the optimum nuclear parton distribution 
functions were determined by the \ 2 analysis of 
DIS and Drell-Yan data. They could be used for 
calculating high-energy nuclear cross sections. 

Second, a possibility of explaining the NuTeV 
sin 2 9w was investigated by the nuclear correc- 
tion difference between u v and d v in the iron nu- 
cleus. Although the contribution to the sin 2 9w 
deviation may not be large at this stage, the dis- 
tribution s v (x) should be investigated by future 
experiments. 

Third, a possible HERMES- type effect was pro- 
posed in the medium and large x regions due to 
the nucleon Fermi motion and binding. Espe- 
cially, we found that the perpendicular nucleon 
motion to the virtual photon direction gives rise 
to the admixture of longitudinal and transverse 
structure functions in the nucleon. Such an effect 
should be tested by electron and neutrino scat- 
tering experiments at large x with small Q 2 . 
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